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ABSTRACT 


The  nonlinear  mechanical  properties  of  a  unidirectional,  glass  fiher- 
epoxy  composite  at  l6i»°F  are  characterized  by  using  creep  auad  recovery 
tests  together  with  a  constitutive  equation  based  on  a  thermodynamic  theory. 
Our  experimental  effort  covers  studies  on  specimens  with  fiber  orientations 
at  0  =  0°,  30°,  45°,  60°  and  90°  with  respect  to  the  uniaxial  loading 
direction.  The  results  for  0  =  0°specimen  show  that  the  composite  is 
linearly  elastic  for  the  load  range  studied,  but  nonlinear  viscoelastic 
behavior  is  observed  for  the  other  orientations.  Guided  by  the  nonlinear 
constitutive  equation,  the  creep  and  recovery  data  are  plotted  on  double- 
logarithmic  paper,  and  the  material  properties  are  found  by  shifting  the 
data  to  form  a  "master  curve"  for  each  fiber  orientation.  Prediction  of 
master  curves  of  0  =  45°  and  60°  from  master  curves  of  0  =  0°,  90°  and  30° 
is  made.  Four  principal  creep  compliances  are  estimated  by  using  the 
master  curves  and  tensor  transformation  relations.  Finally,  we  use  the 
nonlinear  equation  to  predict  strain  response  due  to  multiple-step  loading 
and  unloading. 

The  Appendix  contains  two  abstracts  of  papers  which  were  completed  and 
published  during  the  period  covered  by  this  report. 
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Section  I 


INTRODUCTION 


Nonlinear  viscoelastic  behavior  of  a  unidirectional,  glass  fiber-epoxy 
composite  under  controlled  loading  and  at  approximately  l6U°F  and  21^ 
relative  humidity  is  studied.  The  nonlinear  constitutive  equation  which 
is  used  here  has  been  derived  from  thermodynamic  theory  [l]*  and  is  similar 
to  the  Boltzmann  superposition  integral  of  linear  viscoelastic  theory. 

Special  cases  of  this  nonlinear  theory  have  been  applied  with  limited 
success  by  other  investigators  to  characterize  monolithic  polymeric 
materials.  For  example,  Leaderman  proposed  and  applied  to  fibers  a  so-called 
modified  superposition  principle  (MSP)  [2],  Recently,  Findley  and  Lai  [3] 
used  a  MSP  to  predict  strain  response  of  unplasticized  PVC  to  discontin- 
uously  applied  tensile  and  shearing  stresses.  In  a  more  recent  paper  [It] 
the  various  nonlinear  theories  are  reviewed,  and  the  type  of  theory  used 
here  is  applied  to  experimental  data  on  different  materials  under  controlled 
loading  and  controlled  straining. 

In  Section  II,  linear  and  nonlinear  constitutive  equations  for  uniaxial 
loading  sire  given.  An  equation  for  the  average  octahedral  shear  stress  in 
the  epoxy  matrix  is  also  derived  for  later  reference. 

The  specimen  and  experimental  equipment  are  described  in  Section  III. 
Creep  and  recovery  data  for  various  loads  and  fiber  angles  are  then  given. 

It  is  shown,  using  the  nonlinear  theory,  that  a  graphical  shifting  procedure 


*  Numbers  in  square  brackets  indicate  references  at  end  of  report. 
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can  be  used  to  reduce  both  the  nonlinear  creep  and  the  recovery  data 
plotted  on  double-logarithmic  paper.  This  procedure  provides  master  curves 
which,  for  the  creep  period,  are  simply  the  time-dependent  linear  visco¬ 
elastic  creep  compliances.  Stress-dependent  nonlinear  properties  are 
determined  by  the  amount  of  vertical  and  horizontal  shifting  needed  to  form 
the  master  curves;  we  find  these  properties  depend  primarily  on  the  average 
octahedral  shear  stress  in  the  epoxy.  How  well  data  form  master  creep  and 
recovery  curves  provides  an  indication  of  the  validity  of  the  theory ,  The 
results  reported  herein  do  indeed  verify  the  theory  in  that  good  master 
curves  are  obtained.  Prediction  of  linear  viscoelastic  compliances  at 
0  =  1^5^  and  60^  from  the  curves  of  0  =  0^,  90^  30^,  are  then  made  by 

using  tensor  transformation  relations. 

In  Section  IV  the  principal  creep  compliances  are  estimated.  Finally, 
strain  response  due  to  multiple-step  loading  is  predicted  in  Section  V 
using  the  nonlinear  constitutive  equation  and  the  material  properties  found 
from  creep  and  recovery  data.  Good  agreement  between  theory  and  experiment 
is  achieved  during  most  of  the  time  range  covered. 
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Section  II 


CONSTITUTIVE  THEORY 

l‘.  LINEAR  CONSTITUTIVE  EQUATION 

The  Boltzmann  superposition  integral  can  be  used  to  represent  the 
stress-strain  relation  for  linear  viscoelastic  materials.  Strain  response, 
e,  to  an  arbitrary  stress  input,  o,  is 

e  =  /^  A(t-T)^  dT  (la) 

or ,  equivalently , 


e  =  A(0)a  +  /^  AA(t-T)~  di  (lb) 

o 

where  A(t)  is  the  creep  compliance;  it  is  defined  as  the  strain  output  due 
to  a  unit  step  stress  input,  a  =  H(t),  and 


H(t) 


0  T  <  0 

1  T  >  0 


(2) 


The  range  of  integration  in  (l)  includes  the  stress  discontinuity  that 
occurs  at  t  =  0  in  a  creep  test.  A(0)  is  the  initial  value  of  creep  compli¬ 
ance  and  AA(t)  is  the  transient  component  of  creep  compliance,  AA(t)  = 

A(t)  -  A(0). 

The  linear  viscoelastic  compliance  A(t)  is  a  function  of  the  angle 
between  fiber  and  loading  axis,  9.  Later,  we  shall  use  the  notation  Ag 
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to  bring  out  this  angular  dependence, 

V/hen  a  constant  stress  is  applied  at  t  =  0  and  removed  at  t  = 

viz.  , 


a  =  a^[K(T)  -  H(T-t^)]  (3) 

ve  have  a  so-called  ^'creep  and  recovery  test,” 

Substituting  (3)  into  (lb),  the  stress-strain  relation  for  the  creep  period, 
0  <  t  <  t^jis 

E  =  A(0)a  +  AA(t)o  =  A(t)a  (ii) 

O  GO 

while  during  recovery,  t  >  t^,  it  is 

=  £(t)  -  £(t-t^)  (5a) 

or,  equivalently, 

e  =  [AA(t)  -  AA(t-t_)]o  (5b) 

r  1  o 

In  (5),  is  the  strain  measixred  after  the  removal  of  the  stress,  and  £.(t) 
is  the  creep  strain  for  ^  ^  which  would  exist  if  the  stress  had  not  been 
removed  at  t  =  t^.  The  familiar  relation  between  creep  and  recovery  pre¬ 
dicted  by  equation  (5)  for  a  J.inear  viscoelastic  material  is  shown  in 
Figure  1, 


k 


F'igure  1. 

Relation  Between  Creep  and  Recovery 
for  a  Linear  Viscoelastic  Material 
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2.  NONLINEAR  CONSTITUTIVE  EQUATION 

The  nonlinear  constitutive  equation  derived  in  [l]  is 

t 

e  =  g  A(0)a  +  g,  /  HiAi'p-ii'  )■  ■  ■  ■■-  dT  (6) 

o  o 

where  A(iJ;)  =  A.(o)  +  AA(tp)  is  still  the  linear  creep  compliance  and  ijj 
is  the  so-called  "reduced-time"  defined  "by 

<(j  =  /^  dr'  /a  la(T' )]  (Ta) 

o  ° 

i|<'  =  tpir)  =  dr'/a  [0(1')]  (7b) 

o  ^ 

and  the  material  properties  g^,  g^,  and  are  functions  of  stress  in¬ 
put.  When  the  material  is  loaded  within  the  linear  range ,  ~  ~  ^2 

a^  =  1,  and  equation  (6)  "becomes  the  Boltzmann  superposition  integral  (l). 

The  modified  superposition  principal  proposed  by  Leaderman  [2]  can  be 
derived  from  equation  (6)  by  setting  =  1;  it  is 

e  =  A(0)a  +  AA(t-T)-r^  dx  (8) 

o  dr 

The  thermodynamic  significance  of  changes  in  and  g^  is  that 

they  reflect  third  and  higher  order  stress  dependence  of  the  Gibb’s  free 
energy,  while  arises  from  similar  strong  stress  effects  in  both  entropy 
production  and  free  energy. 

When  a  step  stress  input  a  =  a^[H(T)  -  H(T-t^)]  is  applied,  equation 
(6)  yields, 
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e  =  g  A(0)a  +  g  g  A{~)a  ,  0  <  t  <  t  (9a) 

O  O  i.  3.  o  -L 

Q 

=  g2[AA(!;;)  -  AA(^-iJ^^)  ]0^  ,  <  t  (9b) 


where  in  (9'b) 


Notice  that  e  s  =  g^A(0)a^  +  g^g2AA(iJ^^)a^  immediately  before  the 
load  is  removed,  and  e  =  =  g^AA(\lj^)a^  immediately  after  the  load  is 

removed.  Therefore  the  jump  in  strain  at  t  =  t^  is 

Ae.  =  e,-e'  =  g  A(0)a  +  (g  -l)g  AA(i(»  )a  (lO) 

1  ll  o  o  ^  xo 

It  is  important  to  notice  that  the  strain  Jump  at  t  =  0,  e(0)  =  g^A(0)a^ 
is  not  neces3eu:’ily  equal  to  the  jirnip  at  t  =  t^,  Ae^.  It  is  clear  from 
equation  (lO)  that  these  two  jumps  are  equal  only  when  =  1;  the  linear 
viscoelastic  material  has  the  same  jump  in  strain  during  both  loading  and 
unloading  since  g^  =  1  in  the  linear  range. 

3.  OCTAHEDRAL  SHEAR  STRESS 

The  linear  viscoelastic  creep  compliances  are  the  only  time-dependent 
material  properties  in  the  theory  used  here.  Other  material  properties  are 
stress-dependent  functions,  and  from  our  experimental  results  we  find  that 
they  can  be  expressed  approximately  as  a  function  of  a  single  invariant: 
the  average  octahedral  shear  stress  in  the  epoxy  matrix.  That  the  octahedral 
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shear  stress  is  the  main  function  needed  to  characterize  nonlinear  behavior 
is  suggested  by  raultiaxial  creep  data  on  plastics  and  metals  [^4].  For 
later  reference  purposes,  ve  shall  develop  this  invariant  here  in  terms 
of  applied  stress  and  fiber  angle. 

Octahedral  shear  stress  is  defined  as  the  shear  stress  which  acts 
on  a  plane  whose  normal  makes  equal  acute  angles  with  the  positive  princi¬ 
pal  stress  directions.  It  is  [5] 


where  and  are  the  principal  stresses. 

Consider  now  a  composite  specimen  with  fiber  orientation  6,  under  a 
uniaxial  stress  a,  which  is  shown  in  Figure  2.  The  average  (applied) 
stresses  normal  to  the  fiber,  and  along  the  fiber,  a^,  and  also  the 

average  shear  stress  along  the  fiber,  t^,  are  found  using  the  usual  tensor 
transformation  relations: 

0^  =  o  sin^O  (l2a) 

0^  =  0  cos^0  (l^b) 

Tf  =  ^  sin20  (l2c ) 


—  2L  - 

SPECIMEN 


Figiu*e  2. 

Diagran  for  Derivation  of  Octahedral  Shear  Stress 
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The  average  stress  state  of  the  matrix  element  is  plane  stress;  viz. 

0  =  T  =  T  =0.  Furthermore,  if  ve  assiime  the  fibers  are  so  relatively 

z  xz  yz 

stiff  that  average  stress  on  the  matrix  element  due  to  a 

becomes , 

0=0  (13a) 

y  ^ 


0  =  V  a 

X  e  y 


(13b) 


(13c) 


where  v  is  defined  as  the  ratio  a  it  is  equal  to  the  matrix 

e  XI 

Poisson ratio  in  the  linear  range  only.  However,  it  turns  out  the  result 
is  quite  insensitive  to  and  therefore  the  linear  elastic  value  will  be 
used  in  reducing  our  data. 

Octahedral  shear  stress  is  found  from  the  stresses  (l3)  and  equation 
(11)  to  be 


ilhh) 


We  find  y  has  a  maximum  value  of  h  when  v  =0.5* 

e 

Because  the  fibers  are  not  infinitely  rigid,  is  not  actually  zero. 

Instead,  it  equals  A_.a*  +  A-^a«,and  corrected  stresses  on  the  matrix 
^  11  f  12  f 

element  are 
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(15a) 


=  (E  A,,cot2e  +  E  A,.  +  v)c.  = 


(15b) 


^xy 


(15c) 


where  v  „  =  E  A^TCot^e  +  E  A^-  +  v 
ef  e  11  e  12  e 

is  the  modulus  of  the  epoxy,  and  A^^  and  A^^  principal  composite 
compliance  components  along  the  fiber  and  transverse  to  the  fiber,  re¬ 
spectively. 

The  same  octahedral  shear  stress  form  (l4)  is  obtained  by  using  the 

stress  distribution  (15)  and  equation  (ll),  except  the  effective  Poisson’s 

ratio  V  -  is  used  instead  of  the  epoxy  Poisson's  ratio  v  . 
ef  ^ 

If  we  use  typical  epoxy  properties 


V  =  0.35 
e 


E  =  0.5  X  10°  psi 
e 


and  the  composite  properties  which  we  will  find  later  from  our  experimental 
data  (Aj^^  =  .184  x  10“®  psi"^  and  A^^  =  -.057  x  lO"^  psi“M  it  is  found  that 

Y  =  3.88  when  is  used.  For  the  specimens  studied,  the  biggest  error  in 

Y  due  to  the  use  of  v  instead  of  v  .  is  2.1^  which  occiirs  when  0  =  45°. 

This  small  amount  of  error  in  y  affects  the  octahedral  shear  stress  only 
slightly,  so  that  the  equation  (lU)  will  be  used.  We  should  add  that  a 
quantity 


0  =  T  . 

e  oct 


iJx  = 


+  YT^]' 
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will  actually  be  used  in  reducing  data,  for  convenience.  This  invariant 

is  sometimes  called  the  effective  stress.  It  is  observed  that  a  =  a 

e  f 

when  the  matrix  itself  is  under  a  uniaxial  stress  equal  to  o^. 
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Section  III 


EXPERIMENTAL  WORK  AND  DATA  REDUCTION 


1.  SPECIMEN  AND  EQUIPMENT 

The  contracting  agency  supplied  the  unidirectional »  glass  fiber-epoxy 
composite  material.  It  was  fabricated  by  contacting  eight  plies  for  five 
minutes  and  then  applying  a  pressure  of  25  psi  for  35  minutes.  A  tempera¬ 
ture  of  330^F  was  maintained  throughout  the  pressure  cycle. 

The  composite  has  a  density  of  1.67  g/c.c.,  62.4  weight  percent  fiber, 
and  47.6  volume  percent  fiber. 

Five  sets  of  specimens,  with  6  equal  0,  30,  45,  60,  and  90  degrees, 
were  cut  to  the  nominal  dimensions  of  1/2”  x  l/lO”  x  6”  from  the  same  large 
plate.  Three  specimens  were  cut  at  each  angle.  One  of  them  served  as  a 
dummy  to  compensate  the  thermal  expansion  of  the  strain  gage  used,  and  the 
other  two  specimens  were  both  tested  in  the  test  program  of  one  hour  creep 
and  two  hoiirs  recovery.  Some  difference  was  found  between  the  specimens, 
mainly  in  the  initial  (elastic)  response;  since  the  largest  difference  in 
elastic  response  was  found  to  be  only  3.4^,  which  occured  for  6  =  30^,  data 
from  one  specimen  at  each  angle  was  chosen  for  data  reduction  purposes. 

The  equipment  used  is  shown  in  Figure  3.  A  Develco  Environmental 
chamber  is  shown  inside  the  ARA*  loading  frame  toward  the  right  center. 

Two  Budd  Model  TC-22  Digital  Strain  indicators  are  used  to  record  strains, 
and  they  are  shown  in  the  left  center.  A  Bristol  temperature  and  humidity 
controller  is  on  the  right  side  while  on  the  left  side  is  the  control  panel 

*  Allied  Research  Associates,  Inc. 
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I  till  rd  f  1 


of  the  ARA  Precision  Pnematic  Testing  Machine,  A  nitrogen  tank  is  shown 
at  the  left  corner  to  supply  pneiimatic  pressure  for  the  ARA  Testing  Machine, 
which  has  a  loading  capacity  of  l800  pounds. 

In  Figure  4,  a  close-up  view  of  the  specimen  inside  the  chamber  is 
shown;  the  specimen  is  arranged  vertically  between  grips.  (A  plate  in  the 
background  makes  the  specimen  appear  to  be  bent). 

Two  BLH*  strain  gages  (FAER-25R-35-S13)  are  mounted  on  the  specimen. 
They  are  placed  along  the  longitudinal  direction  of  the  specimen,  one  on 
each  side  and  connected  in  series  to  average  out  the  bending.  For  0  and  90 
degree  specimens,  an  additional  gage  is  placed  along  the  transverse  direc¬ 
tion  to  measure  the  transverse  strain  directly.  Each  gage  has  high  resis¬ 
tance  of  350  Ohms.  Eastman  910  is  used  as  the  bonding  agent,  which  has  a 
maximum  allowable  temperature  of  250^F. 

It  was  found  that  during  the  early  testing  period  of  each  specimen 
the  strain  due  to  a  given  load  decreased  with  successive  tests.  After 
five  to  ten  cycles  repeating  the  loading  and  unloading  process,  the  speci¬ 
men  became  stable  in  creep  compliance.  For  example,  a  stress  of  2690  psi 
was  applied  to  condition  our  60  degree  specimen.  We  found  that  total 
strain  output  under  a  one  hour  loading  creep  test  is  reduced  by  3%  after 
several  cycles. 

This  reduction  in  compliance  may  be  due  to  the  arrest  of  cracks  which 
form  in  the  matrix  under  initial  load  application.  That  cracks  do  indeed 
develop  in  the  rigid  matrix  has  been  discussed  by  Halpin  [6], 


*  Baldwin-Lima-Hamilton  Corporation 
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Figure  4 

Specimen  Inside  the  Environmental  Chamber 
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2.  EQUATIONS  FOR  LINEAR  VISCOELASTIC  COMPLIANCES 

For  a.  linGar  viscoela.s'tic  inELGrial  under  constant  stresses,  and 
in  a  state  of  plane  stress  in  the  x-y  plane,  the  strain-stress  relations 
for  any  fiber  orientation,  6,  relative  to  the  x  axis  are 


(17b) 


Y  —  At-  cf  At— o  A 
'xy  6l  X  62  y  66  xy 


(iTc) 


(in  all  of  our  discussions  the  fibers  are  assumed  to  be  in  the  x-y  jilane.) 
The  elements  AI  are  the  creep  compliances  at  this  particular  value 

^t) 

of  0.  They  form  a  fourth  order  tensor,  so  that  the  relation  between  the 
creep  compliances  of  different  fiber  orientations  6  follows  the  usual 
tensor  transformation  law  [?]. 

For  a  state  of  plane  stress,  only  four  principal  creep  compliances 

A  ,  A  ,  A  and  A^^  are  needed  to  characterize  the  composite;  the  prin- 
11  *  12  22  66 

cipal  directions  are  parallel  to  and  normal  to  the  fibers.  V/hen  the  x-y 
axes  are  these  principle  axes,  the  stress-strain  relations  become 


'x  '  *ll°x  * 


‘y  *2i°x  *  ^aa^y 


(I8b) 


^xv  ^66^> 


(I8c) 


IT 


Pagano  and  Halpin  [8]  derived  an  axial  compliance  of  a  tensile  speci 
men  which  is  clamped  at  both  ends: 


6A’  At.  +  A^\2  ^  -  6A1. 

*0  ^  ^  '  li  J  (19) 

^  °  rAtA>+A>2L_ 

^^11^66  ^11  b2 

where  components  of  the  creep  compliance  tensor 

at  this  angle  0,  while  2L  and  2h  are  the  length  and  width  of  the  specimen, 

respectively;  the  complexity  of  this  equation  results  from  specimen  bending 

when  0  4  0^  ^  4  90^.  If  the  specimen  is  sufficiently  long  compared  to 

its  width  (—  <<  l),  equation  (19)  becomes 
b 

*6  ^  f  “  ''ll 

This  approximation  will  be  used  in  all  of  our  data  reduction  of  linear  and 
nonlinear  data,  and  will  be  verified  numerically  later  in  the  paper. 


3.  0  DEGREE  SPECIMEN 


Negligible  creep  has  been  found  when  loaded  along  the  fibers  (0  =  0^) 
at  164^F  and  21^  relative  humidity,  with  a  stress  range  from  1376  psi  to 
6881  psi.  The  axial  strain  output  due  to  uniaxial  loading  a  =  is 
linearly  elastic,  which  is  shown  in  Figure  5*  The  directly  measured  axial 
compliance,  A^,  is  equal  to  the  principal  creep  compliance 


A^  =  A  =  — 
11  o  a 


=  .18^4  X  10  ^  (psi  M 


The  lateral  strain  ^  due  to  the  axial  stress  a  is  also  found  to 

y 

be  linear  and  independent  of  time;  from  equation  (l8b)  the  princijml  creep 
compliance  A^^  is 

e 

Api  =  ~  .060  X  10  ^  (psi  M 

X 
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The  Poisson’s  ratio  of  the  composite  transverse  to  the  fiber  is 


e 

V  =  -  = 

21  £ 

o 


u.  OFF-ANGLE  SPECIMENS 


The  analytical  expression  for  the  nonlinear  creep  compliance  at  any 


angle  is  obtained  from  equation  (9a), 


"  a  =  ^  SlSyAy-^)  (2l) 

where,  for  constant  stress,  =  t/a 

a 

In  another  paper  [U]  a  method  for  evaluating  all  property  functions 
in  (21)  from  creex)  and  recovery  data  is  given.  Depending  on  whether  or 
not  AA^(i|;)  is  a  power  law,  somewhat  different  approaches  are  needed. 
Notice  that  when  it  is  a  power  law,  viz. 

AAg(i>)  =  cy  (22) 

where  c  and  n  are  constants,  equation  (2l)  yields 


This  shows  that  creep  data  will  provide  information  for  evaluating  the 
combination  ^1^2  turns  out  that 

the  power  law  can  indeed  be  used  to  fit  our  data.  This  power  law  repre¬ 
sentation  has  been  reported  for  many  plastics  a,nd  metals  [9]. 


a.  Initial  Compliance 


In  principle,  creep  data  can  be  used  to  evaluate  g  and  A(0)  in 

o 

equation  (23);  recall  that  g  =  1  under  sufficiently  small  stresses,  and  g 

o  •  o 

therefore  represents  nonlinearity  that  may  exist  in  the  initial  compliance. 
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Howev'6r ,  we  were  unable  "to  measure  strain  at  a  short  enough  time  to  obtain 
the  initial  compliance  directly  from  the  data  in  the  nonlinear  creep  ran^e 
of  behavior.  As  a  result,  the  following  method  was  used  to  estimate  the 
initial  compliance. 

For  each  fiber  angle  studied,  we  found  our  first  strain  reading  (at 

approximately  t  =  3  seconds)  is  linear  in  stress  when  the  stress  is 

sufficiently  small;  this  is  approximately  the  same  range  for  which  recovery 

equation  (5)  is  found  to  hold.  The  slope  of  this  3-second  curve  is  assumed 

to  be  the  initial  creep  compliance  A(0).  Motivated  by  the  nonlinear  creep 

behavior  reported  in  [lO]  for  a  fiber-reinforced  phenolic  resin,  in  that 

g  =1  for  all  applied  stresses,  we  assume  this  is  also  the  case  here. 

^’o 

Graphs  of  the  initial  strain,  e(0)  =  A(0)o  are  shown  in  Figure  j.  (ihe 
solid  points  are  some  measured  vedues,  while  the  open  points  are  esti¬ 
mated  values  using  the  constant  initial  compliance  and  each  stress  value 
applied  in  the  test  program).  The  initial  creep  compliances,  as  given  by 
the  slopes  of  these  lines,  are 

A^q{0)  =  O.UlO  X  10“^  (psi'M 

Aj^^(O)  =  0.600  X  lO”^  (psi  M 

AgQ(O)  =  0.668  X  lO”^  (psi  ^) 

A9o(0)  =  0.6l4  X  10  ®  (psi  M 

It  is  found  that  when  the  initial  compliance  determined  above  is  sub¬ 
tracted  from  the  total  compliance, A^^,  the  power  law  shown  in  (23)  results. 
Even  if  our  estimations  of  initial  compliances  are  not  valid  for  very  short 
times,  they  should  be  acceptable  as  long  as  prediction  of  very  short  time 
response  is  not  required  for  the  engineering  application  of  interest. 

When  0  =  90°,  the  initial  creep  compliance  is  equal  to  the  principal 
compliance  A^g ( 0 ) » 
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^22^^^  ”  ^90^^^  ~  O.61U  X  10~^  (psi  M 


Further,  the  lateral  strain  measurement  e  due  to  a  sufficiently  small 

X 

axial  stress  a  normal  to  the  fibers  (0  =  90^)  gives  the  principal  creep 
compliance  ^sing  equation  (l8a).  For  the  entire  stress  range 

studied,  found  to  be  essentially  independent  of  time  and  stress, 

with  the  value 

£ 

A-,,  =  ~  =  -  0.05^  X  10"^  (psi"M 
12  a  ^ 

The  initial  Poisson* s  ratio  of  the  composite  for  loading  normal  to  fibers 
is  found  experimentally  as 

e 

v,^(0)  =  -  —  =  0.088 
jl2  e 

o 

Theory  shows  that  A^^  =  A^^  [ll].  However,  we  reported  previously 
that  A^^  =  -  .060  X  10  ^  (psi  the  approximately  10^  difference  between 

^12  ^21  attributed  to  experimental  error,  considering  the  very 

small  strains  involved.  The  average  value  will  be  used  in  characterizing 
the  material, 

bp  "  bi  "  ■  (psi'M 

The  initial  compliances  are  plotted  against  fiber  angle  in  Figure  6; 
also  shown  is  the  one-hour  compliance  which  will  be  discussed  later. 

In  this  Figure,  compliances  at  0  =  0^,  90^,  and  30^  were  used  in  fourth- 
order  tensor  transformation  relations  to  predict  compliances  at  G  =  ^5^ 
and  60°.  Also  shown  are  data  for  a  second  set  of  specimens,  which  brings 
out  the  specimen-to-specimen  variation  that  exists,  even  though  they  are 
all  cut  out  of  the  same  plate.  It  is  seen  that  noticeably  better  prediction 
of  the  and  60^  compliances  would  have  been  obtained  if  the  larger 

compliance  at  30^  (solid  point)  had  been  used. 
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Some  of  the  error  in  these  predictions  is  probably  due  to  inhomo¬ 
geneity  in  the  plate,  as  vrell  as  our  method  of  estimating  initial  compliance. 

b.  Reduction  of  Creep  Data 

Most  nonlinearity  and  creep  occurred  for  the  30^  specimen.  Vie  shall 
therefore  use  this  fiber  angle  to  discuss  the  details  of  data  reduction. 

Creep  and  recovery  data  on  this  specimen  are  shown  in  Figure  7.  For 
comparison  purposes  v/e  note  that  the  ultimate  tensile  strength  at  30^  is 
approximately  10300  psi,  and  the  material  is  nonlinear  for  o  >  2000  psi. 

Equation  (21)  can  be  rewritten  as 

and 

logAA  =  log  g  g  +  logM  (t/a  )  (25) 

un  J.  u  a 

where  is  the  transient  component  of  the  nonlinear  creep  compliance, 

AAq  is  the  transient  component  of  the  linear  creep  compliance,  and  we  have 
assumed  =  1. 

We  have  plotted  against  t  on  double-logarithmic  paper 

in  Figure  8.  Nearly  parallel  straight  lines  are  shown  for  different  stress 
levels,  with  only  a  small  slope  difference  existing  at  the  lov/est  stress 
level  (2069  psi).  Linear  viscoelastic  relation  (5)  is  found  to  be  valid 
at  this  lowest  stress,  for  which  the  material  is  considered  to  be  approxi¬ 
mately  linear;  we  shall  assume  the  material  is  linear  for  all  stresses 
0  ^  2069  psi.  Recording  error  for  the  small  strain  output  at  low  stress 
and  the  method  used  to  estimate  initial  compliance  might  have  caused  the 
observed  difference  in  slope. 


2h 
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Figure  7.  Creep— keco very  Curves  =  30' 
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It  is  clear  that  the  data  at  different  stress  levels  can  he  shifted 
along  the  two  axes  to  a  parallel  straight  line  near  the  line  at  the  lowest 
stress  level,  which  is  shown  as  the  broken  line  in  Figure  8.  It  is  the 
so-called  master  curve,  and  is  assumed  to  be  the  transient  component  oi 
the  linear  viscoelastic  creep  compliance. 

The  amount  of  horizontal  (t)  shift  and  vertical  (aA^^)  shift  eq^uals 
log  a  and  log  respectively.  Since  the  lines  at  different  stress 

levels  are  parallel,  straight  lines,  it  is  obvious  that  we  cannot  obtain 
unique  values  for  a  and  shift  of  the  creep  data;  had  curva¬ 

ture  been  found,  the  vertical  and  horizontal  shifts  would  have  been  deter¬ 
mined  uniquely. 

Additional  information  will  be  provided  by  recovery  data  in  order 
to  evaluate  unique  values  of  a^  and  ^^^2^  also  and  g^  separately, 

this  will  be  discussed  in  the  next  subsection. 

Smaller  nonlinearity  is  shown  in  Figure  9  for  0  =  90^.  Parallel, 
straight  lines  are  also  found  for  the  creep  compliance  data  at  different 
stress  levels  for  0  =  90  ,  as  well  as  0  =  ^5  and  60  . 

c.  Reduction  of  Recovery  Data 

The  nearly  parallel,  straight  lines  on  double-logarithmic  paper  shown 
in  Figures  8  and  9  indicate  that  the  power-law  form  can  represent  the 
linear  creep  compliance;  namely, 

Aq  =  Aq(0)  +  AAg(i|;)  =  Ag(0)  +  (26) 


27 


where  tj;  =  is  stress-reduced  time  and  c  and  n  are  independent  of  stress 

and  time. 

Substituting  this  power  law  (26)  into  equation  (9b)  yields 


and  we  derive  the  equation  for  normalized  recovery  strain 


— [(1  +  a  X)”  -  (a  A)^] 
Ae^  a  a 


X  =  -- 


is  the  transient  component  of  strain  existing  immediately  before  the 
stress  is  removed. 

Rewriting  equation  (28)  we  have. 


log(-^|-)  =  -log  +  log[(l  +  -  (a^^)^]  (29) 

It  is  clear  from  equation  (29)  that  a  graphical  shifting  procedure 

can  be  applied  to  recovery  data  to  determine  and  for  different  stress 

levels.  For  0  =  30°,  it  is  estimated  that  n  =  .145  from  Fijrpire  3,  A 

e 

reference  curve  of  normalized  recovery  strain  “ —  is  plotted  against  A 

A^i 

on  double-logarithmic  paper  by  letting  =  g^^  =  1  in  equation  (29).  This 
reference  curve  is  the  normalized  recovery  strain  for  a  linear  viscoelastic 
material;  it  is  shown  as  the  solid  line  in  Figure  10.  Four  normalized 
recovery  curves  of  different  stress  levels  are  also  shown  in  broken  lines. 

It  appears  that  they  can  be  s)iperposed  on  the  reference  (linear  viscoelastic) 
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curve  by  translating  them  along  two  axes,  the  amount  of  horizontal  (a) 

e 

shift  (to  the  left)  and  vertical  (•^“)  shift  (up)  is  equal  to  -  log  a^ 
and  log  g^,  respectively. 

Since  the  curves  shown  in  Figure  10  are  not  straight  lines,  unique 
values  of  a^  and  gj^  can  be  obtained  by  the  shifting  technique.  After 
shifting,  some  of  the  data  points  for  four  different  normalized  recovery 
curves  are  shown  as  solid  points  on  the  reference  curve.  The  ability  to 
shift  the  data  to  the  reference  curve  provides  a  check  on  the  theory. 

The  normalized  recovery  curve  for  the  lowest  stress  level  (2069  psi) 
is  not  shown  in  Figure  10  since  it  was  found  to  lie  approximately  on  the 
curve  of  3h60  psi,  and  not  on  the  reference  curve.  However,  this  does 
not  mean  the  material  is  nonlinear  at  2069  psi.  Recall  that  the  actual 
slope  of  the  transient  component  of  creep  compliance  at  2069  psi  is  not 
the  same  as  the  master  curve  in  Figure  8;  by  using  the  actual  slope  of 
this  creep  compliance  (n  =  0.18?)  for  predicting  the  reference  recovery 
curve,  we  found  that  the  measured  recovery  curve  at  2069  psi  lies  approxi¬ 
mately  on  this  new  reference  curve;  this  result  shows,  therefore,  that 
the  material  is  essentially  linear  at  the  lowest  stress  of  2069  psi. 

Given  the  values  of  log  a^  from  the  recovery  Figure  10,  the  creep 
curves  in  Figure  8  are  shifted  this  distance  to  the  right .  The  additional 

amount  of  vertical  shift  needed  to  place  each  creep  curve  on  the  master 
curve  (dashed  line)  gives  the  value  of  log  Knowing  g^  from  I'igure  10, 

we  then  determine  the  value  of  at  each  stress. 

This  shifting  procedure  has  been  carried  out  for  each  stress  level 
and  fiber  angle  tested.  The  nonlinear  properties  are  plotted  in  Figures 
11,  12  and  13  against  the  "effective  stress,"  equation  (l6),  which  is 
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Fiber  Angie 


Figure  12. 


O  30 


proportional  to  the  average  octahedral  shear  stress  in  the  matrix.  The 

constant  value  v  =  0.35  'w'as  used  to  evaluate  y  =  3.38,  as  discussed  in 
e 

Section  II~3.  The  fact  that  most  of  the  data  fall  on  a  single  curve  in 

each  Figure  indicates  that  is  the  principal  factor  causing  nonlinearity. 

Since  >  1  in  Figure  13 >  the  jximp  of  strain  when  t  =  t^  is  larger 

than  that  when  t  =  0;  i.e.  >  ^^^(0)0^  in  equation  (lO). 

The  linear  dependence  of  log  a  on  a  seen  in  Figure  11  at  high 

o  e 

stresses  is  especially  interesting.  This  is  the  same  behavior  observed 
earlier  for  a  glass-reinforced  phenolic  [U];  the  fibers  were  in  an  orthog¬ 
onal  array,  so  that  shear  stress  along  the  fibers,  rather  than  was 
found  to  be  the  principal  invariant.  This  linearity  in  Figure  11  for  the 
larger  values  of  implies  that 

.  -Ba 

a  »  Ae  e 
a 

where  A  and  B  are  positive  constants.  Exponential  dependence  is  suggested 
by  Eyring's  rate  process  theory  for  high  stresses  when  all  nonlinearity  is 
due  to  a^[9],  but  this  molecular  model  is  too  simple  to  account  for  and 

^2  * 

The  transient  components  of  recovery  data  are  plotted  against  linear 
time  scales  in  Figures  1^  and  15  for  the  largest  stresses  applied  at 
6  =  30^  and  90^.  The  accuracy  of  the  recovery  prediction  is  an  indication 
of  how  close  the  shifted  recovery  data  in  Figure  10  are  to  the  solid  line. 
The  modified  superposition  theory  prediction,  equation  (8),  was  also  made; 
note  that  equation  (5a)  can  be  used  instead  to  make  this  prediction  since 
it  is  valid  for  nonlinear  behavior  (g^  4  T)  when  =  g^  =  =  1.  Total 

strain  using  (5a)  is  shown  in  Figure  7.  These  results  indicate  that  the 
modified  superposition  theory  of  Leaderman  predicts  recovery  which  is  too 
rapid. 
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Figure  \U*  Net  Creep-Recovery 
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Figure  IS.  Het  Creep-Recovery  Curve  (6*=  9C°) 


Cl*  Prediction  of  Creep  Compliances 

For  the  specimens  0  =  90^^,  ^5^,  and  6o^,  we  found  that  the  material 
is  approximately  linearly  viscoelastic  up  to  stresses  of  2079  psi,  2069  psi 
and  13^5  psi,  respectively.  For  the  linear  range,  fourth  order  tensor 
transformation  relations  will  be  used  to  predict  compliances  at  6  =  4$^ 
and  60*^  given  data  at  0  =  0^,  90^,  and  30^.  In  order  to  predict  compliances 
at  higher  stresses,  one  would  simply  bring  in  the  values  of  g^,  and 
along  with  (9a);  if  all  data  points  in  Figures  11  and  12  were  on  a  single 
curve,  rather  than  being  somewhat  scattered,  the  predictions  in  the  non¬ 
linear  range  would  be  as  accurate  as  in  the  linear  range. 

First,  we  record  the  linear  creep  compliances  found  experimentally. 

They  are 

A  =  .184  X  lO"*^  (psi"M 

o 

=  [.6ll+  +  X  10~^  (psi'U 

=  [.1+10  +  X  10*6  (psi"e 

(30) 

A)^^  =  [,600  +  X  10“^  (psi'M 

AgQ  =  [.668  +  .lOTti^'^"^^]  X  lO'S  (psi“M 

where  if;  is  stress-reduced  time  in  hours;  recall  for  the  linear  range  =  1, 
so  that  ij;  =  t. 

The  transient  components  of  the  creep  compliances  (AA  =AiiJ;)  -  aIo)  = 

U  o  6 

cip^)  have  been  plotted  against  if;  on  double-logarithmic  paper  in  Figure  l6. 

Since  A^^  are  fourth  order  tensor  components,  the  standard  transforma¬ 
tion  relations  [?]  give 

Aq  =  (m^  -  3m^n^)A^  +  (n^  -  1/3  ^  (3l) 

where 

m  =  COS0,  n  =-.sin0 
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Figure  16.  'transient  Component  of  Linear  Viscoelastic  Creep  Compliances 


The  transient  parts  of  the  predicted  and  A’  are  shoTO  in  broken 

oo 

lines  in  Figure  l6  to  compare  with  the  experimental  results  -  compare  also 
equations  (33)  with  (30).  For  further  comparison,  1  hour  data  of  Figure  l6 
are  plotted  against  fiber  angle,  G,  in  Figure  6. 

Part  of  the  observed  error  is  probably  due  to  experimental  errors  and 
the  lack  of  material  homogeneity  of  the  sheet  from  which  the  specimens  were 
cut;  also,  small  errors  in  the  original  estimates  of  initial  compliances 
produce  relatively  large  errors  in  the  transient  components,  and  could  he 
one  cause  of  the  differences  in  slopes  in  Figure  l6. 
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Section  IV 


PRINCIPAL  CREEP  COMPLIANCES 


It  was  mentioned  previously  in  Section  III-2  that  the  principal  creep 
compliances  A^^^  and  A^^^  are  directly  determined  frora  axial  strain  measure¬ 
ments  under  uniaxial  step  stress  input  along  the  fibers  and  normal  to  the 
fibers,  respectively;  and  the  third  principal  creep  compliance,  can  be 

measured  by  a  strain  gage  placed  transverse  to  the  loading  direction  on  a 
0^  or  90*^  specimen.  They  are  found  to  be  independent  of  stress  and  have 
the  values, 

=  .184  X  10~^  (psi"M 

^22  ^  bo  '  ^  X  10“6  (psi~M  (3,^) 

b2  "  bl  (psi”M 


The  fourth  principal  creep  compliance  A^g  can  be  calculated  by  using 
the  data  of  A-,..,  A^^. ,  A  and  A  ,  and  the  standard  tensor  transformation 

11  dtl  12  oO 

relation  [7]: 


‘66  =  r  -  “u  -  2^12  -  I  '‘22 


(39) 


The  resulting  four  principal  creep  compliances  are  plotted  against 
log  ip  in  Figure  17.  These  four  compliances  completely  characterize  linear 
viscoelastic  behavior  under  all  states  of  plane  stress. 

Earlier  we  said  equation  (20)  is  only  an  approximation  because  of  speci¬ 
men  bending  that  exists  with  off-angle  specimens.  Now  we  want  to  check  this 
approximation.  By  using  the  principal  data  and  tensor  transformation 


hi 


rinci p  al 


which  would  he  measured  if 


relations  we  can  calculate  the  compliance 
the  material  really  had  the  principal  creep  compliances  above.  It  appears 
that  the  approximation  is  very  good,  with  the  directly  measured  value 
as  shown  in  Table  1. 

Table  I 

Comparison  of  the  measured  and  calculated  A^q. 

ij;  Measured  x  10®  Calculated  A^q  x  10® 

(Mrs)  (psi)"^  (psi)  ^ 


0 

.410 

.407 

1 

.472 

.468 

20 

.507 

.502 

200 

.547 

.541 

^3 


Section  V 


.^4ULTIPLE-STEP  LOADING 


The  nonlinear  theory  has  been  successfully  applied  to  creep  and 
recovery  data  reduction  and  prediction  in  the  previous  Sections.  This 
theory  will  now  be  used  to  predict  response  under  inult iple-step  loadinf:^ 
and  then  compared  to  test  data  for  the  6  =  30^  specimen. 

The  constitutive  equation  (6)  is  used  to  make  this  prediction, 


e  fr^  A(0)a  +  g-  /  AA(D-4;  ^  (ix 

o  In  dT 


ip  =  f  dt/a 


The  two-step  loading  and  two-step  unloading  program  shown  in  Figure  l8 


was  chosen. 


0  q  t 


Figure  l8 

Multiple-Step  Loading  Program 


The  stress  input  function  has  the  form, 

■  gpp(T)  +  (g^o^  -  g2^j^)H(T-tj^) 


-(RoCt,)  -  g^a  )H(T-t  )  -  glo  H(T-t_) 


The  related  nonlinear  properties  a^,  and  have  a  similar  form;  for 
example , 

a  =  1  +  (a^  -  1)H(t)  +  (a^  -  aMH(T-t  ) 

O  0  O  O  X 

-(a2  -  aMH(T-t„)  -  (a^  -  l)H(T-t_)  (37) 

0  0^0  j 

where  the  superscript  on  a  material  property  indicates  the  particular 
stress  at  which  it  is  to  be  evaluated. 

Substituting  the  stress  history  (36)  into  equation  (6),  and  using  the 
related  nonlinear  properties,  e.g.  (37),  yields, 
for  0  <  t  <  t^ 

e  =  g^A(0)a^  +  AA(tJj)a^  (38a) 

ip  =  t/a^ 

for  t^  <  t  <  tg 

£  =  g2A(0)a2  +  g^[g2^^AA(iii)  +  (g|a^.,  -  g^a^)  AA('j;-ii'^)  ] 
ip  =  +  (t-t^)/a^ 

=  (t-t^)/a^  (38b) 

for  tg  <  t  <  t^ 

e  =  glA(o)0^  +  g^[g^o^AA(<(;)  +  (g^a^  -  gJa^)AA(4'-'l'^) 

-(ggOg  -  g^a^)AA(i|;-i;<2)]  (38c) 

f  =  t^/al  +  (t2-tp)/a2  +  (t-t2)/al 

<P-ip^  =  (t-t2)/aj 
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for  <  t , 

e  =  g^a^AA(iJ;)  + 

-(g^Og  -  gla^)M(;>-ij^^)  -  g^a^AA(i;)-i{<^)  (38d) 

i|;  =  t^/al  +  (t2-t^)/a2  +  (t^-t2)/al  +  (t-t^) 

=  (t2-t^)/a2  +  (t^-t^)^^  +  (t-t^) 

4i-<p^  =  {t^-t^)/a}^  +  (t-t^) 

The  actual  stress  and  time  values  used  are  =  3^^60  psi,  =  69^0  psi, 

and  t^  =  1/2  hr,  t^  1  hr,  t^  =  1.9  Irir.  The  test  was  made  twice  on  the 

same  30^  specimen  at  l6^4^F  and  21a  relative  huraidity  used  to  characterize 

the  material.  The  linear  creep  compliance  and  the  stress-dependent  nonlinear 

properties,  which  we  obtained  previously  from  data  reduction  procedures,  are 

used  in  equations  (30)  to  predict  the  strain  response.  The  results  for  the 

second  repetition  of  this  load  history, only  slightly  different  from  those 

of  the  first  load  history,  are  shown  in  Figure  19  and  seen  to  agree  very 

well  with  the  experimental  observations  except  for  the  third  region 

(t,^  <  t  <  t  ).  This  figure  shows  that  the  agreement  is  much  better  than 
2  3 

achieved  using  linear  theory. 

The  data  in  Figure  19  were  obtained  several  weeks  after  the  specimen 
was  characterized.  Therefore,  some  of  the  error  can  probably  be  attributed 
to  the  long  rest  period  and  resulting  change  in  properties. 


k6 


Fiber  Angle  6  =  30® 
Temperature  *  164  ±  2®F 
Rel. Humidity  *2I±I% 


d  d 

(%)3‘Nivais 
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Figure  19.  Strain  Response  to  Multiple-Step  Loading 


Section  VI 


CONCLUSIONS 


An  appreciable  amount  of  nonlinearity  in  the  transient  component  of 
strain  was  found  to  exist  in  all  specimens,  except  for  those  under  loading 
parallel  to  the  fibers.  The  nonlinear  constitutive  equation  (6)  derived 
from  thermodynamic  theory  [l]  has  been  used  successfully  to  characterize 
the  glass  fiber-epoxy  composite  under  uniaxial  loading;  although  only  a 
uniaxial  equation  has  been  given  here,  multiaxial  relations  have  been 
derived  [l].  Linear  viscoelastic  compliances  and  nonlinear  properties 
were  evaluated  by  applying  graphical  shifting  procedures  to  creep  and  re¬ 
covery  data.  All  three  nonlinear  pronerties,  ,  g,..  and  a  were  found  to 

’  -^1’  a 

be  a  function  of  essentially  only  one  mechanical  invariant:  the  average 
octahedral  shear  stress  in  the  epoxy  matrix. 

Such  a  result  seems  to  indicate  uniaxial  creep  and  recovery  tests 
may  be  sufficient  to  characterize  the  composite  for  multiaxial  stress  appli¬ 
cations.  However,  the  conclusion  is  very  preliminary  at  this  time,  as 
other  functions  in  [l]  may  enter  and  confirmation  from  a  multiaxial-stress 
test  program  is  needed. 

The  nonlinear  theory  used  in  this  report  is  applied  more  easily 
than  the  familiar  multiple- integral  representation  in  characterizing  a  non¬ 
linear  viscoelastic  solid.  The  single-integral  form  makes  it  relatively 
easy  to  use  in  engineering  stress  analysis  applications  as  well.  Moreover, 
the  only  time-dependent  functions  which  appear  in  the  constitutive  theory 
are  the  linear  viscoelastic  creep  compliances. 


Some  error  was  observed  in  predicting  the  angular  dependence  of  the 
compliances  and  strain  response  to  multiple-step  loading.  It  is  believed 
that  possible  sources  of  the  error,  in  addition  to  spec imen-to- specimen 
property  variations  and  experimental  error,  are  the  method  used  to  estimate 
initial  compliance  and  a  gradual  change  in  properties  observed  for  a  single 
specimen  under  repeated  loading  and  unloading. 

Motivated  by  this  finding,  a  new  data  reduction  procedure  has  been 
developed  without  restricting  =  1  [12];  recovery  data  is  utilized  in 
estimating  the  initial  compliance.  This  new  method  provides  better  agree¬ 
ment  with  respect  to  angular  dependence  of  the  compliance  but  does  not 
appreciably  change  the  response  under  multiple-step  loading. 

The  effect  of  temperature  on  material  properties  is  not  covered  in 
the  present  report,  but  [12]  does  include  some  results  at  various  tem¬ 
peratures  . 
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APPENDIX  I 


On  the  Application  of  a  Thermodynamic  Constitutive 
Equation  to  Various  Nonlinear  Materials* 

R.  A,  Schapery 

ABSTRACT 

The  author’s  thermodynamic  constitutive  theory  for  nonlinear  visco¬ 
elastic  behavior  is  extended  to  account  for  rate-independent  plastic  flow 
and  for  nonlinear  creejj  with  strong  stress-dependence.  The  resulting 
constitutive  equation  is  then  applied  to  several  different  materials  under 
small  and  large  strains  in  order  to  bring  out  the  significance  of  the  non¬ 
linear  functions,  and  to  obtain  a  check  on  the  theory.  History  effects 
are  accounted  for  by  means  of  a  single  integral  which  is  very  similar  to 
the  Boltzmann  form  in  linear  theory,  and  therefore  permits  easy  application 
to  the  various  loading  and  straining  conditions  studied.  On  the  basis  of 
physical  arguments  and  some  example  applications  to  glassy  and  polycrystalline 
materials,  including  fiber-reinforced  plastic,  it  is  suggested  that  nonlinear, 
multiaxial  response  of  these  solids  can  often  be  described  by  linear  visco¬ 
elastic  equations  with  an  intrinsic  time-scale  that  is  affected  by  mechani¬ 
cal  invariants  and  temperature;  for  example,  the  classical  three-dimensional 
stress-strain  equations  of  rate-independent,  incremental  plasticity  and 
nonlinear  creep  of  metals  are  shown  to  be  included  in  this  description. 

Only  a  slightly  more  involved  form  is  needed  to  characterize  a  noiaber  of 
unfilled  and  filled  soft  polymers  under  finite  strain. 

*  Presented  at  the  lUTAM  Symposium  on  Thermoinelasticity,  East  Kilbride, 
Scotland,  June  ?6-28,  I968.  (To  be  published  in  Proceedings.)  Work 
was  originally  published  as  Purdue  University  Report  No.  AA&ES  68-^ 

(June  1968). 
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APPFJNDIX  II 


Thermal  Expansion  Coefficients  of  Composite 
Materials  Based  on  Energy  Principles* 

R.  A.  Schapery 

ABSTRACT 

Bounds  on  effective  thermal  expansion  coefficients  of  isotropic  and 
anisotropic  composite  materials  consisting  of  isotropic  phases  are  derived 
by  emi^loying  extremum  principles  of  thermoelasticity.  Inequalities  between 
certain  approximate  and  exact  forms  of  the  potential  and  complementary 
energy  functionals  are  first  established.  These  inequalities  are  then 
used  in  conjunction  with  a  new  method  for  minimizing  the  difference  between 
upper  and  lower  bounds  in  order  to  derive  volumetric  and  linear  thermal 
expansion  coefficients.  Application  is  made  to  two-  and  three-phase  iso¬ 
tropic  composites  and  a  fiber-reinforced  material.  It  is  found  for  some 
important  cases  that  the  solutions  are  exact,  and  take  a  very  simple  form. 
V/e  also  show  conditions  under  wiiich  the  rule-of-mixtures  and  Turner's 
equation  can  be  used  for  thermal  expansion  coefficients.  Finally,  simple 
methods  for  extending  all  results  to  viscoelastic  composites  are  indicated. 


*Published  in  J.  Composite  Materials,  Vol.  2,  pp.  3'60-h0h  (July  I968). 
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